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Abstract 

We model the underdoped cuprates using fermions moving in a background with local antifer- 
romagnetic order. The antiferromagnetic order fluctuates in orientation, but not in magnitude, 
so that there is no long-range antiferromagnetism, but a 'topological' order survives. The normal 
state is described as a fractionalized Fermi liquid (FL*), with electron-like quasiparticles coupled 
to the fractionalized excitations of the fluctuating antiferromagnet. The electronic quasiparticles 
reside near pocket Fermi surfaces enclosing total area x (the dopant density), centered away from 
the magnetic Brillouin zone boundary. The violation of the conventional Luttinger theorem is 
linked to a 'species doubling' of these quasiparticles. We describe phenomenological theories of the 
pairing of these quasiparticles, and show that a large class of mean-field theories generically dis- 
plays a nodal-anti- nodal 'dichotomy': the interplay of local antiferromagnetism and pairing leads 
to a small gap near the nodes of the d-wave pairing along the Brillouin zone diagonals, and a large 
gap in the anti-nodal region. 



I. INTRODUCTION 



The nature of the ground state in the underdoped regime of the hole-doped cuprate 
superconductors remains a central open issue. Angle resolved photoemision spectroscopy 
(ARPES) and scanning tunneling microscopy (STM) have been the main tools to explore 
such a regime. In both probes, an unexpected angular dependence of the electron spectral 
gap function has been revealed: a 'dichotomy' between the nodal and anti-nodal regions of 
the Brillouin zone in the superconducting stated, Specifically, this dichotomy is realized by 
deviations in the angular dependence of the gap from that of a short-range d-wave pairing 
amplitude ~ (cos k x — cos k y ). 

This paper will describe the superconducting instabilities of a recently developed modeP 
of the normal state of the underdoped cuprates based upon a theory of fluctuating local 
antiferromagnetic order^33. A related normal state model of fluctuating antiferromagnets 
has been discussed by Khodas and TsvelikP^, who obtained results on the influence of spin- 
wave fluctuations about the ordered state similar to ours 6 . These results have been found 
to agree well with ARPES observations^^. Another approach using fluctuating antifer- 
romagnetism to model the underdoped cuprates has been discussed recently by Sedrakyan 
and Chubukov^. We will also connect with the scenario emerging from recent dynamical 
mean- field theory (DMFT) studies™ 3 . 

The theory of Ref. El describes the normal state in the underdoped regime as a fractional- 
ized Fermi liquid (FFL or FL*), although this identification was not explicitly made in that 
paper. So we begin our discussion by describing the the structure of the FL* phase. 

The FL* phase is most naturally constructed^^ using a Kondo lattice model describing 
a band of conduction electrons coupled to lattice of localized spins arising from a half-filled 
d (or /) band. The key characteristics of the FL* are {%) a 'small' Fermi surface whose 
volume is determined by the density of conduction electrons alone, and (m) the presence 
of gauge and fractionalized neutral spinon excitations of a spin liquid. In the simplest 
picture, the FL* can be viewed in terms of two nearly decoupled components, a small Fermi 
surface of conduction electrons and a spin liquid of the half-filled d band. The FL* should 
be contrasted from the conventional Fermi liquid, in which there is a 'large' Fermi surface 
whose volume counts both the conduction and d electrons: such a heavy Fermi liquid phase 
has been observed in many 'heavy fermion' rare-earth intermetallics. Recent experiments 
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on YbRh 2 (Sio.95Geo.o5)2 have presented evidence^ for an unconventional phase, which could 
possibly be a FL*. 

A concept related to the FL* is that of a "orbital-selective Mott transition"^ (OSMT), 
as discussed in the review by Vojta 24 . For latter, we begin with a multi-band model, like 
the lattice Anderson model of conduction and d electrons, and have a Mott transition to 
an insulating state on only a subset of the bands (such as the d band in the Anderson 
model). The OSMT has been described so far using dynamical mean field theory (DMFT), 
which has an over-simplified treatment of the Mott insulator. In finite dimensions, any such 
Mott insulator must not break lattice symmetries which increase the size of a unit cell, for 
otherwise the state reached by the OSMT is indistinguishable from a conventionally ordered 
state. Thus the Mott insulator must be realized as a fractionalized spin liquid with collective 
gauge excitations; such gauge excitations are not present in the DMFT treatment. With a 
Mott insulating spin liquid, the phase reached by the OSMT becomes a FL*. 

Returning our discussion to the cuprates, there is strong ARPES evidence for only a single 
band of electrons, with a conventional Luttinger volume of 1+x holes at optimal doping and 
higher (here x is the density of holes doped into the half-filled insulator). Consequently, the 
idea of an OSMT does not seem directly applicable. However, Ferrero et al™ argued that 
an OSMT could occur in momentum space within the context of a single-band model. They 
separated the Brillouin zone into the 'nodal' and 'anti-nodal' regions, and represented the 
physics using a 2-site DMFT solution. Then in the underdoped region, the anti-nodal region 
underwent a Mott transition into an insulator, while the nodal regions remained metallic. 
A similar transition was seen by Sordi et al. in studies with a 4-site cluster^. While these 
works offers useful hints on the structure of the intermediate energy physics, ultimately the 
DMFT method does not allow full characterization of the different low energy quasiparticles 
or the nature of any collective gauge excitations. 

We turn then to the work of Ref. El who considered a single band model of a fluctuating 
antiferromagnet. Their results amount to a demonstration that a FL* state can be con- 
structed also in a single band model, and this FL* state will form the basis of the analysis of 
the present paper. The basic idea is that the large Fermi surface is broken apart into pockets 
by local antiferromagnetic Neel order. We allow quantum fluctuations in the orientations of 
the Neel order so that there is no global, long-range Neel order. However, spacetime 'hedge- 
hog' defects in the Neel order are suppressed, so that a spin liquid with bosonic spinons and 
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a U(l) gauge-boson excitation is realizecP^^. Alternatively, the Neel order could develop 
spiral spin correlations, and suppressing Z 2 vortices in the spiral order realizes a Z 2 spin 
liquid with bosonic spinons^^. The Fermi pockets also fractionalize in this process, and we 
are left with Fermi pockets of spinless fermions; the resulting phase was called the algebraic 
charge liquicP^ (ACL). Depending upon the nature of the gauge excitations of the spin 
liquid, the ACL can have different varieties: the U(1)-ACL and SU(2)-ACL were described 
in RefsJ-^, and Z 2 -ACL descends from these by a Higgs transition involving a scalar with 
U(l) charge 2, as in the insulatoi * 27 * 28 l 

Although these ACLs are potentially stable phases of matter, they are generically suscep- 
tible to transformation into FL* phases. As was already noted in Ref. El there is a strong 
tendency for the spinless fermions to found bound states with the bosonic spinons, leading 
to pocket Fermi surfaces of quasiparticles of spin S — 1/2 and charge ±e. Also, as we 
will review below, there is a 'species-doubling' of these bound states^^, and this is crucial 
in issues related to the Luttinger theorem, and to our description of the superconducting 
state in the present paper. When the binding of spinless fermions to spinons is carried to 
completion, so that Fermi surfaces of spinless fermions has been completely depleted, we are 
left with Fermi pockets of electron/hole-like quasiparticles which enclose a total volume of 
precisely x holes^. The resulting phase then has all the key characteristics of the FL* noted 
above, and so we identify it here as a FL*. The U(1)-ACL and Z 2 -ACL above lead to the 
conducting U(1)-FL* and Z 2 -FL* states respectively. Ref. El presented a phenomenological 
Hamiltonian to describe the band structure of these FL* phases. Thus this is an explicit 
route to the appearance of an OSMT in a single-band, doped antiferromagnet: it is the local 
antiferromagnetic order which differentiates regions of the Brillouin zone, and then drives 
a Mott transition into a spin liquid state, leaving behind Fermi pockets of holes/electrons 
with a total volume of x holes. 

We should note here that the U(1)-FL* state with a U(l) spin liquid is ultimately unstable 
to the appearance of valence bond solid (VBS) order at long scales®'. However the Z 2 -FL* is 
expected to describe a stable quantum ground state. The analysis of the fermion spectrum 
below remains the same for the two cases. 

Phases closely related to the U(1)-FL* and Z 2 -FL* appeared already in the work of 
Ref. This paper examined 'quantum disordered' phases of the Shraiman-Siggia modeP^, 
and found states with small Fermi pockets, but no long-range antiferromagnetic order. The 
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antiferromagnetic correlations where either collinear or spiral, corresponding to the U(l) 
and Z 2 cases. However, the topological order in the sector with neutral spinful excitations 
was not recognized in this work: these spin excitations were described in terms of a 0(3) 
vector, rather than SU(2) spinor description we shall use here. Indeed, the topological order 
is required in such phases, and is closely linked to the deviation from the traditional volume 
of the Fermi surfaces) 20 * 21 ! 

We also note another approach to the description of a FL* state in a single band model, 
in the work of Ribeiro, Wen, and RarP^^I They obtain a small Fermi surface of electron- 
like "dopons" moving in spin-liquid background. However, unlike our approach with gapped 
bosonic spinons (and associated connections with magnetically ordered phases), their spinons 
are fermionic and have gapless Dirac excitation spectra centered at (±7r/2, ±7r/2). 

We will take the U(1)-FL* or Z2-FL* state with bosonic spinon spin liquid as our model 
for the underdoped cuprates in the present paper. We will investigate its pairing properties 
using a simple phenomenological model of (f-wave pairing. Our strategy will be to use the 
simplest possible model with nearest-neighbor pairing with a <i-wave structure, constrained 
by the requirement that the full square lattice translational symmetry and spin-rotation 
symmetry be preserved. Even within this simple context, we will find that our mean-field 
theories of the FL* state allows us to easily obtain the 'dichotomy' in the pairing amplitude 
over a very broad range of parameters. We also note that the pocket Fermi surfaces of 
the FL* state will exhibit quantum oscillations in an applied magnetic field with a Zeeman 
splitting of free spins, and this may be relevant to recent observations^. 

We mention here our previous workP^El on pairing in the parent ACL phase. These 
papers considered pairing of spinless fermions, while the spin sector was fully gapped: this 
therefore led to an exotic superconductor in which the Bogoliubov quasiparticles did not 
carry spin. In contrast, our analysis here will be on the pairing instability of the FL* state, 
where we assume that the fermions have already bound into electron-like quasiparticles, as 
discussed above and in more detail in Ref. M The resulting Bogoliubov quasiparticles then 
have the conventional quantum numbers. 

Our primary results are illustrated in Fig. [T] We also show a comparison to a conven- 
tional state with co-existing spin density wave (SDW) and d-wave pairing, and to recent 
experiments. The left panels illustrate Fermi surface structures in the normal state. The 
right panels show the angular dependence of the electron gap in the superconducting states: 
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(a) (b) 




FIG. 1: Color online: Our new results for the FL* phase (second row), compared with the Hartree- 
Fock/BCS theory (top row) and experiments (bottom row). The left panels illustrate Fermi surface 
structures in the normal state. The right panels shows the angular dependence of the electron gap 
in the superconducting states, (a) Spectral weight of the electron in the normal state with SDW 
order at wavevector K = (tt,tt). Here we simply apply a potential which oscillates at (vr,7r) 
to the large Fermi surface in the overdoped region, (b) Minimum electron gap as a function of 
azimuthal angle in the Brillouin zone. The full (red) line is the result with a pairing amplitude 
~ (cosA; x . — cos k y ) co-existing with SDW order, while the dashed (black) line is the normal SDW 
state, (c) Spectral weight of the electron in the FL* state, with parameters as in Fig. |3j note that 
the pocket is no longer centered at (7r/2,7r/2). (d) Spectral gap functions in the superconducting 
(full (red) line) and normal (dashed (black) line) states of Fig. [3j (e) The Fermi pocket from a 
ARPES experiment)^ related observations appear in Refs. fT2)13l (f) The dichotomy of the spectral 
gap function from the observations of Ref. [3j See the text for more details. 
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for each angle 6, we determine the minimum electron spectral gap along that direction in 
the Brillouin zone, and plot the result as a function of 9. 

The results of the traditional Hartree-Fock/BCS theory on SDW order and <i-wave pairing 
appear in (a) and (b). The SDW order has wavevector K = (it, tt), and the c?-wave pairing is 
the conventional (cos k x — cos k y ) form. In the normal state, the Fermi pocket is centered at 
the magnetic Brillouin zone boundary, as shown in (a). An important feature of this simple 
theory is that the state with co-existing SDW and <i-wave pairing has its maximum gap at 
an intermediate angle, as shown in (b): this reflects the "hot spots" which are points on the 
Fermi surface linked by the SDW ordering wavevector K. No experiment has yet seen such 
a maximum at an intermediate angle. 

One set of our typical results for the FL* theory are shown in (c) and (d). As it was shown 
in the previous worlP, the normal state in (c) shows a Fermi pocket which is clearly not 
centered the magnetic zone boundary (at (±7r/2, ±7r/2)); furthermore, its spectral weight 
is not the same along the Fermi surface, and has a arc-like character. At the same time, 
(d) shows the angular dependence of the electron gap in the superconducting state; unlike 
the SDW theory, this FL* state has a pairing amplitude which is a monotonic function of 
angle and has its maximum at the antinodal point. It also shows the "dichotomy" in the 
gap amplitude between the nodal and anti-nodal regions. For the purpose of comparison, we 
illustrate two experimental results in (e) and (f). Clearly, our mean-field theory can provide 
reasonable explanation for the experimental data in both the normal and superconducting 
states, and we believe it is a candidate for the under-doped cuprate materials. 

The structure of this paper is following. In Section [TTJ we introduce the normal state 
Hamiltonian for the fermions, and investigate the symmetry transformations of possible 
pairings. We classify possible pairings which preserve full square lattice symmetry, and 



introduce a low energy effective pairing Hamiltonian. In Section III, spectral gap functions 
for various cases are illustrated assuming d x 2_ y 2 wave pairing. It is shown that our model 
can reproduce the dichotomy behavior, and we compare our theory with the YRZ model 
proposed by Yang, Rice and ZhangJM^ 21 , anc [ the re lated analyses by Wen and Le d 44 1 45 1 . For 
completeness, it is shown that U(l) gauge fluctuation can mediate the needed d wave pairing 
in Appendix O 
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II. EFFECTIVE HAMILTON! AN 



The basic setup of the FL* state has been reviewed in some detail in Refs. IqITiTH and so 
we will be very brief here. The starting pointP ^ 31 * 46 ^ is to transform from the underlying 
electrons Ci a to a rotating reference frame determined by a matrix R acting on spinless 
fermions ip p . 

Cia = Kpi^p- (!) 

R ap is a SU(2) matrix with a =f, I for spin index, p = ± for gauge index, and we parame- 
terize 

R = p ~f] (2) 

with \z,j\ 2 = 1 . In the ACL state, the bosonic z a and the fermionic ip p are assumed to be 
the independent quasiparticle excitations carrying spin and charge respectively. Then we 
examined the formation of bound states between these excitations. A key result was that 
was a "doubling" of electron-like quasiparticles, with the availability of two gauge neutral 
combinations, 

F ia ~z ia ip i+ , G ia ~ e a pz*pij)i-.. (3) 

This doubling is a reflection of the 'topological order' in the underlying U(l) or Z2 spin 
liquid; it would not be present e.g. in a SU(2) spin liquid^. The F ia and the G{ a will be 
the key actors in our theory of the FL* phase here. Their effective Hamiltonian is strongly 
constrained by their non-trivial transformations under the space group of the Hamiltonian, 
which are listed in Table [H 

From these symmetry transformations, we can write down the following effective 
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TABLE I: Transformations of the lattice fields under square lattice symmetry operations. T x : 
translation by one lattice spacing along the x direction; Rty% ' ^0° rotation about a dual lattice 
site on the plaquette center (x — > y, y — > —x); I£ ual : reflection about the dual lattice y axis 
(x — > —x,y —7- y); T: time-reversal, defined as a symmetry (similar to parity) of the imaginary 
time path integral. Note that such a T operation is not anti-linearP 

HamiltoniarP 

Htot — H o + Hi n t 

Ho = -^^^io + ^^ + A^(-l)^(ilF fa -GiG fa ) 

ij i 

(FlG ja + GlF ja ). (4) 

i<j 

Here tij is taken to be similar to the bare electron dispersion, characterizing the Fermi surface 
in the over-doped region; A is a potential due to the local antiferromagnetic order; and 
is the analog of the Shraimain-Siggia termP^ which couples the two species of electron-like 
quasiparticles F and G to each other; it is this term which is responsible for shifting the 
center of the pocket Fermi surfaces in the normal state away from the magnetic Brillouin zone 
boundary. Hi nt is the invariant interaction Hamiltonian: there could be many interaction 
channels, which induces superconductivity of the (F, G) particles, such as negative contact 
interaction, interaction with other order parameters, and the gauge field fluctuation. In this 
paper, we do not specify particular interaction and we assume that pairings are induced. 
Then we focus on studying properties of possible pairings and their consequences on physical 
quantities such as spectral gaps. In Appendix [Cj we illustrate one possible channel to achieve 
such superconductivity. 
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For some of our computations, it is more convenient to use an alternative basis for the 
fermion operators 

The C and D fermions have the same space-group transformation properties as the physical 
electrons. Then, the Hamiltonian becoms 

..^H'f*' * ] M (6 ) 

We chose energy dispersion's forms following the previous workP, with e(k) a Fourier trans- 
form of tij and e(k) a Fourier transform of t^, and K = (tt, tt): 

e(k) = — 2ti(cos k x + cos k y ) + 8t 2 cos k x cos k y — 2t 3 (cos2k x + cos 2k y ) 

e(k) = — 1 — 2ti(cos k x + cos k y ) + 8^2 cos k x cos k y — 2t 3 (cos 2k x + cos2A;j / ) 

e c (k) = e(k) + e(k) - n 

e d (k) = e(k + K) - e(k + K) — //. (7) 

The C and D particles have spin and electric charges like electrons. Therefore, any linear 
combination can be a candidate for the physical electron degree of freedom. In the previous 
workP, we matched the C particles to the electrons of large Fermi surface state without 
antiferromagnetism; following this, for simplicity we will take the C to be the physical 
electron, but our results do not change substantially with other linear combinations. Then 
the D particles are emergent fermion induced by fluctuating SDW order. Note that the 
C, D particles live in the full first Brillouin zone of the square lattice, and not the magnetic 
Brillouin zone. 

Issues related to the Luttinger theorem were discussed in previous workP^E]. The total 
area of the Fermi pockets described by Hq is precisely x, the dopant hole density. Here the 
area is to be computed over the full first Brillouin zone of the square lattice, as the full 
square lattice symmetry is preserved by our model. Also note that our phenomenological 
Hamiltonian H has been designed to apply only to low energy excitations near the Fermi 
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surface. However, rather than focusing on these momentum space regions alone, considera- 
tions of symmetry are far simpler if we define the dispersion in real space on the underlying 
square lattice, as we have done here. For this somewhat artificial lattice model, as discussed 
in Ref. El the total fermion density on each site % is 

E ( C L C ^ + D U D J) = E ( F L F ^ + G{ a G ha ) = 2-x (8) 

a a 

The traditional Luttinger theorem measures electron number modulo 2, and so it should 
now be clear that occupying the independent electron states of the lattice H will yield a 
Fermi surface with the desired area of x. 

Before proceeding further, let us review the above discussion. We started our theory with 
electrons in one band, and considered spin density wave fluctuation. The strong fluctuation 
induced particle fractionalization, and bound states whose degree of freedoms are doubled 
appeared. The resulting phase is nothing but the FL* we introduced above. Therefore, the 
ACL phase provides a natural way to connect the FL* with one band theory. 

To study superconductivity of the FL* phase, let us consider invariant pairing operators. 
With the (F, G) particles, there are many possible combinations in principle. However, it is 
more convenient to work in terms of the C and D particles because they transform just like 
electrons under the symmetry operation. So we can write down the 4 pairing operators 

O c A (iJ) = e^C^Cw, O d A (iJ)=e^D i>a D j , p 

Of{i,j) = s^C i>a D jjfi , 0%(i,j)=e a PD i , a CM (9) 

Note that we only consider even parity pairing, and there are only three pairings, 
O c , O d , O cd + O dc (see Appendix g. 

III. SPECTRAL GAP 

Throughout this paper, we assume that all pairings are d wave, more specifically, d x 2_ y 2. 
The assumption of the d wave pairings can be realized by the gauge fluctuation (see Ap- 
pendix [Cj) or by other channels like conventional spin density wave fluctuations. Then, with 
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the pairing amplitudes as in Eq. (J9]), we can write down the mean field Hamiltonian 
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where A c is the Fourier transform of O^, A,i is the Fourier transform of O^, and Ax is the 
Fourier transform of O c £ + O d £. For their wavevector dependence we take the forms 

A c (k) A d (k) A x (k) 



^cO 



A 



d0 



A 



cos k x — cos fc^ 



xo 



(11) 



where A c0 , A d0 and Axo are the respective gap amplitudes. 

In principle, we could determine these pairing amplitudes from solving a set of BCS-like 
self-consistency equations. However, in the absence of detailed knowledge of the pairing 
interactions, we will just treat the A c o, A^o and Axo as free parameters. In other words, 
we are in the deep superconducting phase with adjusted parameters. Then our task is to 
study spectral gap behaviors with given band structures and pairings. More technically, the 
Green's function of the C particle, which determines the electron properties, are studied 
focusing on the pole of the C particles' Green's function. The pole basically contains infor- 
mation about the electron's dispersion relation, and its minimum determines spectral gap 
properties. The latter is defined as the minimum gap along a line from the Brillouin zone 
center at an angle 9: thus the nodal point is at 9 = it/A, and the anti-nodal point at 9 = 0. 

Although we have three free gap parameters, our results are quite insensitive to their 



values. For simplicity we will mainly work (in Sections III A and III B ) with the case with 
a single gap parameter A c0 ^ 0, and others are set to zero A d0 = Axo = 0. We will briefly 



consider the case with multiple gap parameters in Section |III C[ and find no significant 
changes from single gap case. 
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FIG. 2: Color online: Spectral gap functions and the Fermi surfaces with the Case I, (*2 = 0.15*i, 
* 3 = -0.3t 2 , *i = -0.25*!, * 2 = 0, * 3 = 0, to = -0.3ti, At = -0.6*i, A = 0.4*i. (a) The spectral 
weight of the electron Green function with the relaxation time r*i = 200. (b) Fermi surfaces of 
e c (k) (dashed inner (red)) , e^(k) (dashed outer (blue)), and the eigenmodes (thick (black)) of Hq. 
The dotted line is the magnetic zone boundary, (c) The spectral gap function with and without 
A c . The dotted (black) line is for the normal case. The thick (red) line is for superconducting state 
with A c o = 0.1*i. (d) The spectral gap function with and without A^x- The dotted (black) line 
is for the normal case. The thick (green) line is for the superconducting state with Axo = 0.1*i. 
The dashed (blue) line is the superconducting state with A^o = 0.1*i. 

A. Single Gap : case I 

We consider the case with * 2 = 0.15*i, t 3 = —0.3*2, h = — 0.25*i, * 2 = 0, *3 = 0, 
*o = —0.3*i, fJ> = —0.6*i, A = 0.4*i in Fig. [2] In (a), the calculated spectral weight of the C 
particle is illustrated following the previous paperP The shape is obviously pocket-like, but 
its spectral weight depends on position on the Fermi surface. In (b), we illustrate the bare 
energy Fermi surfaces and their eigenmode Fermi surface. Note that the two bare energy 
bands (e C) d(k)) are different from the usual SDW formations with Brillouin zone folding. In 
the latter, there is only one electron band, and SDW onset divides the Brillouin zone two 



13 



pieces (e(k),e(k + K)). But in our case, the two bands have different energy spectrums 
of the electron-like particle (C) and the emergent particle (D). And A determines mixing 
energy scale between the C and D particles. 

In (c), the spectral gap function with and without a given pairing, A c is illustrated. Near 
the node, it is obvious that the pairing gap contributes to the spectral gap in a d wave pairing 
way as expected. However, between the node and anti-node, there is a huge peak. The peak 
position is nothing but the mixing point between C, D particles. Therefore, the peak exists 
whether there is a pairing or not. Near the anti-node, the spectral gap is bigger than the 
near-node's but much smaller than the mixing point peak. It indicates there is tendency 
to make electron pockets near the anti-node. For example, if we decrease the magnitude 
of A, which basically represent the mixing energy scale, then the gap near the anti-node 
becomes smaller, and eventually the electron-like pockets appears near the anti-node with 
the pre-existing hole type pockets. (See the Appendix) Note that this situation is formally 
the same as the pairing with the SDW fluctuation mediating pairing case (see Fig. [I]). 
The "hot spot" between the node and the anti-node has the largest gap magnitude, which 
corresponds to our mixing point. Such a spectral gap behavior is not the experimentally 
observed one. Therefore, we cannot have the needed dichotomy near the anti-node in this 
case; the anti-nodal gap is always smaller than the one of the maximum mixing point. 
Following the similar reasoning, the experimentally observed dichotomy does not appear 
in the conventional SDW theory unless additional consideration beyond mean-field theory 
is included. In (d), we illustrate other pairing cases (A d;X )- As we can see, the role of 
the pairings are similar to the conventional one (A c ), and qualitatively they are the same. 
Therefore, it is not possible to achieve the observed dichotomy by considering the exotic 
pairings. They cannot push the maximum peak of the normal state to the anti-nodal region. 

The message of this calculation is simple. With the band structure we considered here, the 
observed dichotomy in the spectral gap function cannot be obtained, even though the normal 
state can explain experimentally observed Fermi surface structures. Moreover, it also implies 
that it is difficult to explain the observed dichotomy with the Hartree-Fock/BCS mean-field 
theory of the Fermi liquid. 

However, we now show how our FL* theory gets a route to explain the dichotomy below. 
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FIG. 3: Color online: Spectral gap functions and the Fermi surfaces for the Case II (*2 = 0.15*i, 
£ 3 = -0.3£ 2 , h = -0.25ti, t 2 = 0, £3 = 0, £ = -0.3£i, y, = -0.8£i, A = 0.6*i). Note that the 
only change from Fig. [2] is in the values of \x and A. (a) The spectral weight of the electron Green 
function with the relaxation time rt\ = 200. (b) Fermi surfaces of e c (dashed inner (red)) , Ed 
(dashed outer (blue)), and the eigenmodes (thick (black)) of Hq. The dotted line is the magnetic 
zone boundary, (c) The spectral gap function with and without A c . The dotted (black) line is for 
the normal case with A c = 0. The thick (red) line is the superconducting state with A c o = 0.1£i. 
(d) The spectral gap function with and without A^x- The dotted (black) line is for the normal 
case. The thick (green) line has Axo = 0.1£i. The dashed (blue) line has A^o = 0.1£i. 

B. Single Gap : case II 



In Fig. |3|we illustrate the case with £2 = 0.15*i, *3 = —0.3*2, £1 = — 0.25*i, *2 = 0, £3 = 0, 
* = —0.3*i, A* = ~~ 0.8*i ) A — 0.6*i. These parameters are as in Section HI A, except that 



the values of \i and A have changed. As we discuss below, this changes the structure of the 
dispersion of the 'bare' C and D particles in a manner which leaves the normal state Fermi 
surface invariant, but dramatically modifies the spectral gap in the superconducting state. 

As we can see in (a), the calculated spectral weight of the C particle is qualitatively 
the same as the Fig. [2js. The shape is obviously pocket-like, and its spectral weight also 
depends on position of the Fermi surface similarly. Therefore, in the normal state, there is 
no way to distinguish the two cases because the low energy theory are all determined by 
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the Fermi pocket structures. However, in (b), the bare energy Fermi surfaces of e c (k) and 
€d(k) are clearly different from the previous one's. Even though the bare Fermi surfaces look 
unfamiliar, they are irrelevant for the observed Fermi surface which is determined by the 
eigenmodes of Hq (black line), and which is qualitatively the same as the Case I. 

We illustrate our spectral gap behavior with and without the pairing, A c , in (c), which 
was already shown in the introductory section. Without the given pairing, the normal state 
has the finite gap less region where the pockets exist, and there is a stable spectral gap in the 
anti-node. It is easy to check the anti-nodal gap depends on the mixing term, A, between 
the C and D particles. With the pairing, the Fermi pockets are gapped and only the node 
remains gapless. The spectral gap function has expected d wave type gap near the node, and 
the observed dichotomy is clearly shown. Therefore, the origin of the two gaps are manifest; 
the nodal gap is obviously from the C particle pairing and the anti-nodal gap is originated 
from the mixing term, which is inherited from the spin-fermion interaction term. In (d), 
we illustrate other exotic pairings (A^x)- As we can see, role of the pairings are similar to 
the conventional pairing (A c ), and qualitatively they are the same. So, there is no way to 
distinguish what pairings are dominant only by studying spectral gaps. 

Now let us compare our results to the ones of the YRZ modeP^^. In the YRZ model, 
based on a specific spin liquid model, the pseudo-gap behavior is pre-assumed by putting 
an explicit d x 2_ y 2 gap function in the spectrum, which means the characteristic of the anti- 
nodal gap is another input parameter. With the two d wave gaps (pairing and pseudo-gap), 
the experimental results were fitted. 

In our FL* theory, the anti-nodal gap behavior is determined by the interplay between 
A and the bare spectrum e Ci d(k) Indeed, the pseudo-gap corresponding term, A, is s wave 
type in terms of YRZ terminology. The A term represents local antiferromagnetism, and 
this 'competing' order which plays a significant role in the anti-nodal gap. The parameter 
A is just input for making the Fermi pockets in the normal state with other dispersion 
parameters. As mentioned before, it explains the distinct origins of the nodal and anti- 
nodal gaps. Also, although our theory contains other pairings, A^x, we did not need that 
freedom to obtain consistency with experimental observations. 

Of course, non-local terms of A could be considered. And it is easy to show that the 
d x 2_ y 2 like terms are not allowed because of the rotational symmetry breaking. Putting the 
non-local A term is secondary effect, and we do not consider it here. 
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FIG. 4: Multiple gaps. The left panel is the same as the Fig. [2] with two superconducting gaps 
Ado = 0.3ii and A c o = 0.1t±. And the right panel is the same as the Fig. [3] with A^o = 0.3ii and 
A c o = O.lti. In both, the dashed (green) line is with the two gaps. And the plain and dotted lines 
are the same as the previous plots. 

C. Multiple gaps 

So far, we have only considered the cases with one pairing gap. Of course, multiple gaps 
are possible and we illustrate possible two cases in Fig. |4j which contain A c ^ with the two 
normal band structures. Here, we choose the same phase in both pairings. The spectral 
gap behaviors are not self-destructive, which means the magnitude of spectral gap with two 
pairings is bigger than the one with the single pairings. One comment is that even multiple 
gaps do not change qualitative behavior of the spectral gap functions, which means that the 
Case I could not have the observed dichotomy even with the multiple gaps. 

In Fig. [5j two pairings with the opposite sign are illustrated. Clearly, we can see the 
self-destructive pattern with the same gap magnitudes. Even a node appears beyond the 
nodal point. Therefore, it is clear that the relative phase between two pairings plays an 
important role to determine the gap spectrum. 

IV. CONCLUSIONS 

This paper has presented a simple phenomenological model for pairing in the underdoped 
cuprates, starting from the FL* normal state described in Ref. El This is an exotic normal 
state in which the Cu spins are assumed to form a spin liquid, and the dopants then occupy 
states with electron-like quantum numbers. A key feature of this procedure^, is that there 
is a 'doubling' of the electron-like species^ available for the dopants to occupy: this appears 
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FIG. 5: Multiple gaps with the relative phase difference. Details are the same as Fig. [3j The red 
line is for two superconducting gaps with the same sign gaps, A^o = O.lti and A c o = O.lii. And 
the green line is for the opposite sign gaps, A^o = — O.lti and A c o = 0.1t\. 

to be a generic property of such doped FL* states. 

Our previous worlP showed how this model could easily capture the Fermi surface struc- 
ture of the underdoped normal state. In particular, a mixing between the doubled fermion 
F and G species from the analog of the 'Shraiman-Siggia' terrrP^led to Fermi pockets which 
were centered away from the antiferromagnetic Brillouin zone boundary. 

Here we considered the paired electron theory, assuming a generic <i-wave gap pairing of 
the cos k x — cos k y variety. Despite this simple gap structure, we found two distinct types of 
electron spectral gaps in this case, illustrated in Figs. [2] and [3] The distinction arose mainly 
from the strength of a parameter, A, determining the strength of the local antiferromagnetic 
order. 

For weaker local antiferromagnetic order, and with a normal state Fermi surface as in 
Fig. [2^, the angular dependence of the gap had a strong maximum near the intermediate 
"hot spot" on the underlying Fermi surface. A similar structure is seen in the traditional 
Hartree-Fock/BCS theory of SDW and <i-wave pairing on a normal Fermi liquid, and this 
structure is incompatible with existing experiments. 

For stronger local antiferromagnetic order, we were able to maintain the normal state 
Fermi surface as in Fig. [3^,, but then found a gap function which had the form shown in 
Figs. |3]3,d, which displays the 'dichotomy' of recent observations. Thus in this theory, it is 
the fluctuating local antiferromagnetism which controls the dichotomy. 

Finally, we compare our theory with model proposed by Yang, Rice, and Zhang^^, and 
closely related results of Wen and Le d 44 * 45 *. Their phenomenological form of the normal state 
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electron Green's function has qualitative similarities to ours^, but there are key differences 
in detail: 

(i) The 'back end' of the YRZ hole pocket is constrained to be at (tt/2, vt/2), while there is 
no analogous pinning in our case. 

(m) The electron spectral weight vanishes in the YRZ theory at (tt/2, 7r/2), while our theory 
has a small, but non-zero, spectral weight at the back end. 

(Hi) Our theory allows for a state with both electron and hole and pockets, while only hole 
pockets are present in the YRZ theory. 

These differences can be traced to the distinct origins of the 'pseudogap' in the two theories. 
Our pseudogap has connections to local antiferromagnetism which fluctuates in orientation 
while suppressing topological defects. Pairing correlations also play an important role in 
the pseudogap, but these are neglected in our present mean-field description: these were 
examined in our previous fluctuation analyses of the ACLpED The YRZ pseudogap is due 
to a d-wave 'spinon pairing gap' in a resonating valence bond spin liquid. All approaches 
have a similar transition to superconductivity, with a <i-wave pairing gap appearing over 
the normal state spectrum, and a nodal-anti-nodal dichotomy: thus any differences in the 
superconducting state can be traced to those in the normal state. 

The differences between our normal state theory with bosonic spinons, and other work 
based upon fermionic spinons^^ESSSl become more pronounced when we consider a transi- 
tion from the normal state to a state with long-range antiferromagnetic order. In our theory, 
such a transition is naturally realized by condensation of bosonic spinons, with universal 
characteristics discussed earlier^^Sl. Such a natural connection to the antiferromagnetically 
ordered state is not present in the YRZ theory. 
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FIG. 6: Color online: Spectral gap functions and the Fermi surfaces for the Case II (*2 = 0.15*i, 
t 3 = -0.3*2, h = -0.25ii, * 2 = 0, t 3 = 0, to = -0.3ti, M = -0.6*i, A = 0.25*i). Note that the only 
change from Fig. [3] is in the value of A. (a) The spectral weight of the electron Green function with 
the relaxation time rt\ = 200. (b) Fermi surfaces of e c (dashed inner (red)) , £d (dashed outer 
(blue)), and the eigenmodes (thick (black)) of Hq. The dotted line is the magnetic zone boundary, 
(c) The spectral gap function with and without A c . The dotted (black) line is for the normal case 
with A c = 0. The thick (red) line is the superconducting state with A c o = 0.05*i. (d) The spectral 
gap function with and without A^x- The dotted (black) line is for the normal case. The thick 
(green) line has Axo = 0.05*i. The dashed (blue) line has A^o = 0.05*i. 

Appendix A: Electron pockets 



We consider the case with t 2 — 0.15*i, £ 3 = — 0.3t 2 , ^i = — 0.25*i, t 2 = 0, t 3 = 0, 
to = —0.3*1, A* = —0.6*i, A = 0.25*i in Fig. pi These parameters are as in Section HI A 



except that the value of A has lowered. In other words, the 'bare' spectrums are the same, 
but electron pockets near the anti-node appear due to the low mixing term. 

As we can see in (a), the calculated spectral weight of the C particle shows the hole and 
electron pockets with different spectral weights. We illustrate our spectral gap behavior with 
and without the pairing in (c) and (d). Without pairings, the normal state has the finite 
gap less region where the pockets exist, and there is an intermediate region peak similar to 
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FIG. 7: Color online : Spectral gap behaviors varying with A. The thick(red), dotted (black) and 
dashed(green) lines are for X/ti = 0.4,0.3,0.25 with the same pairing magnitude, Ax = 0.05ti. 

the SDW case. With pairings, the Fermi pockets are gapped and only the node remains 
gapless. The spectral gap function shows similar behavior as in our case I. In Fig. [7j spectral 
gap function varying the the mixing term is illustrated to see the evolution of the dip near 
the anti-node. 

We note that electron pockets can also appear in the YRZ formulation, but have very 
different shapes^. 

Appendix B: Invariant pairings 

There are four combinations of invariant pairing terms of the F and G: 

0£(i,j) = e aP (F ija F jtfi + G i>a G^) 

o a A {i,j) = ni i-\y ■ , >(F,.j-- J .>-c,j: J .,) 

0\{i,j) = e^i-iy^iG^F^-F^G^) (Bl) 

In Table [TTJ we illustrate the transformation of various pairing terms. The four pairings 
havean interesting exchange symmetry. Obviously 0^ B have even under the exchange oper- 
ation. If we consider nearest neighbor sites, it is easy to show that O b A is even and 0\ 
is odd under the exchange. Therefore, for the d x i_ y i symmetry, the O a does not contribute 
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Ddual 


rdual 


1 


e aP F a Fp 


e a ^G a Gp 


e al3 G a G/3 


e al3 G a Gj3 


-e^FlFl 




c aB TP TP 


c a0 TP TP 


c a0 TP TP 




e aP F a G p 


e a ^G a F^ 


e al3 G a Fp 


e al3 G a Fp 


-e^G\Fl 


e a ?G a Fp 


e aP F a Gp 


e aP F a G p 


e afs F a Gp 


-e^FlGl 



TABLE II: Transformations of the pairing terms. We suppress the lattice index(i, j) before and after 
transformations. Note that the Time Reversal column (T) contains (— ) term and the conjugate 
partner also have the (— ) sign. 

to pairings. 

The conversion between the two representations are as follows: 



O c A (i,j) 


— s a ^Ci jCe Cj^ — 


\(o£ + o»)(i,j) 




O d A (t,j) 


= e al3 D ita Dj^ = 




Ol)(i,j) 


O c ^,j) 


= E al3 Ci :Ct Dj i p = 


\(0 A + O b A )(i,j) 




O d A c (i,j) 


= e a/3 D iQ Cjfi = 


(_1)Ax+Ay 

2 


O b A )(i,j), 



where Ax + Ay is coordinates' difference between two particles, for example, zero for the s 
wave and one for the d wave. 



Appendix C: Pairing Instability 

In this section, we introduce one way to achieve the d wave instability from the gauge 
fluctuation. There could be many other channels to induce the d wave channel such as 
"conventional" SDW fluctuation, so this section shows possibility of obtaining the desired 
pairings. 

To constrain the Hamiltonian, let us consider symmetry transformations of the field 
strengths associated with the U(l) gauge field of the CP 1 model describing the z a spinons 
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rdual 
J x 
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-B 


+B 


-B 


c 

<^x 


c 

^x 


c 


+£-x 


c 

^x 


Cy 


Cy 


+s x 


Cy 


Cy 




t x $> 









TABLE III: Symmetry transformations of the U(l) field strength of the CP 1 model, and of the 
fermion field $ = (FG) T . 

in Table 




Pauli matrix, r(er), is defined in the (F,G) (spin) space. The only invariant coupling up to 
the second order derivatives is^ 

<S 7 = 7 / £■ &T y {V)^ 

J T,X 

= -il I S-(FlVG a -GiVF a ). (C2) 

J T,X 

It is interesting to note that this coupling is precisely the geometric phase coupling between 
the antiferromagnetic and valence bond solid (VBS) order parameters discussed recently in 
Ref. H21 The electric field is the spatial component of the skyrmion currrent in the Neel 
state, and it couples here to a fermion operator which has the same quantum numbers as the 



spatial gradient of the phase of VBS order; thus Eq. (C2) corresponds to the spatial terms 
in Eq. (3.8) in Ref. [521 Here we see that the electric field couples to a 'dipole moment' in 
the fermions. 

We can also look for a coupling between the magnetic field, B, and the fermions. There 
is no coupling up to the second order derivatives of fermionic fields. The main reason for 
the absence is that rotation and inversion transformations have opposite signs acting on the 
magnetic field. If we go beyond the second order derivative, we can find a coupling to the 
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FIG. 8: Vertex correction of a pairing channel. The wavy line is for the propagator of the electric 
fields, and the plain (dotted) line is for the C(D) particle's. Note that the coupling with the 
electric field (filled dot) contains the momentum component. Here we represent the C particle 
pairing vertex renormalization. 

magnetic field such as 

S B = lB [ W\%-tf)(d x d v )<T*V. (C3) 

J T,X 

This term is also one associated with the geometric phase between the antiferromagnetic 
and VBS orders, corresponds to the temporal term in Eq. (3.8) in Ref. |52j 
The fluctuations of the gauge field are controlled by the action 

S * = ^rE/ ^[^ E (k,e n )\S\ 2 + U B (k,e n )\B\ 2 ], (C4) 

where IIe and U B are polarization functions from the matter fields. Because of the non- 
minimal coupling between the electric and magnetic fields and the fermions, there is no 
screening, and these polarization functions are just constants at low momenta and frequen- 
cies. Also, although the bosonic spinons do couple minimally to the electromagnetic field, 
they are gapped and also yield only a constant contribution to the polarizations. 
With the C, D representation, the coupling term to the electric field becomes 

<S 7 = -7 / £{q, n)-k[Dl(k + q + Q,co + Q)C a {k, to) - tfjk + q,u + Q)D a {k + Q, u)] . 

It is manifest that C and D particles are only mixed with the finite momentum Q difference. 
Let us consider pairing vertex 

V pmrmg = nk)C{,Cl k ^ + S d (k)D^Dl k>i + h.c. (C5) 

k 

To see the superconducting instability, we need to evaluate the vertex correction of the 
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pairing channel such as the diagram in Fig. |8j The presence of the A requires numerical 
evaluations. Instead of considering numerical calculations, let us turn off the mixing term, 
A, and see which pairings are preferred with approximations. We will discuss about the 
non-zero mixing term later. 

The renormalized pairing vertex of C particles is 

s c (k) ren ~ ^)- 7 2 4/( fc + «r / -TT-a> ( C6 ) 

J £UJ u 2 + ej 

where I /He is the constant electric field propagator. As usual, we assume that the inte- 
gration is dominant near Fermi surfaces and the k 2 becomes k 2 F . Also we extract the gap 
function of D particles out of the integration. The factor Md is the density of states of 
D particles. Note that the minus sign in front of the second term is from the momentum 
dependence of the interaction and the relative sign of the gap functions. Likewise, the D 
particle pairing correction is 

5\k + Q)ren ~ S*{k + Q) - >fk* F J°(k)£ [ (C7) 

n £ J euJ u 2 + e 2 c 

In both equations, the last integrals show the usual BCS type logarithmic divergence. We 
can determine momentum dependence of the pairings with these equations. If we assume s 
wave pairings, then the corrections become negative and the renormalized pairings become 
suppressed. On the other hand, d wave pairings can change the sign of the integration and 
enhance the superconductivity. Such momentum dependence results from the momentum 



dependent vertex term in Eq. (C2) with a given relative pairing sign. In the gauge exchange, 
the momentum dependence plays the same role as spin-exchange in the usual d wave BCS 
pairing. 

So far, we have fixed the relative sign between the two pairings by hand. Our calculation 
indicates possibility of d wave pairings, but the channel of the instability can vary with 
changing the relative pairing sign. There could be fully gapped pairing with opposite signs, 
s ± . 

Evaluating the vertex corrections, we have assumed no mixing term, A, at the lowest 
approximation. Now let us turn on the mixing term. Then, the Fermi surfaces of the two 
particles start mixing and details of the Fermi surfaces change. Of course, (C, D) pairings 
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can be mixed by A. But the mixing point is first gapped out and the Fermi surfaces become 
pockets. So there is no significant pairing mixing by A and we can treat pairings separately. 
Details of the Fermi surface change, but we can argue that pairing channels remain intact 
at low energy. 
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